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ABSTRACT 

It is suggested that the abundances of neutral 
non-polar interstellar molecules unobservable by 
radioastronomy can be .systematically determined by radio 
observation of the protonated ions. As an example, 

4. # 

observed N 2 H column densities are analyzed to infer 
molecular nitrogen abundances in dense interstellar 
clouds. The chemistries and expected densities of 
the protonated ions of C> 2 , C 2 , C0 2 , C 2 H 2 anc ^' CIi 4 are 
then discussed. Microwave transition frequencies 

. j, 

for HC0 2 and C 2 H 2 are estimated^ and a preliminary 
astronomical search for HC0 2 *" is described. 
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I . INTRODUCTION 


The discovery of an intricate interstellar 
chemistry has been almost entirely an achievement 
of radio astronomy. Large, high density clouds — 
the regions where molecules exist in appreciable 
abundance — are optically opaque; it is therefore 
important for an understanding of the chemistry of 
the interstellar medium to rely on radio observations 
and to extend these to new classes of molecules. In 
this paper, we shall explore the feasibility of deter- 
mining by indirect means the abundances of non-polar 
neutral species which are not in themselves directly 
observable by radio techniques. We are especially 
interested in species such as N 2 , 0 2 , C 2 , CC> 2 , C 2 H 2 , 
and CH^, all of which may be present to an appreciable 
extent'. in dense interstellar clouds. In addition, we 
are interested in an indirect determination of the H 2 0 
abundance because the amount of water present cannot 
be obtained from the maser transition observed. It is 
necessary to determine the abundances of all of these 
species if the chemistry of dense clouds is every to 
be fully understood. 
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An abundance for any species can be obtained 
theoretically from a model of molecule formation 
and depletion if that model is complete with respect 
to all processes that affect the species of inter- 
est to an appreciable extent. (See for example, 
Herbst and Klemperer 1973, hereafter abbreviated as 
HK) . In this paper, a much less complete analysis 
is utilized, in particular, we examine the abund- 
ances of the protonated polar ions N 2 H , 0 2 H , 

C 2 H + , HOCO* H 3 0 + , C 2 H 3 J "' anc ^ CH 5 + P re 3icted relative 
to their neutral precursors . It is assumed that 
these ionic species are produced and destroyed by 
gas phase ion-molecule and ion-electron reactions. 
The rate coefficients of the relevant reactions are 
frequently known from laboratory measurements with 
high accuracy (10-20%) and in general may be approxi 
mated by simple theoretical techniques (HK) . For 
exothermic ion-molecule reactions, the rate coef- 
ficient is almost always near the Langevin value of 
-9 3 -1 

2 x 10 cm sec whereas for ion-electron disso- 
ciative recombination reactions, the rate coeffici- 
ent is usually ~ 1-2 x 10 cm" 3 sec - (HK) . From 
the appropriate kinetic equations, it is possible 
to deduce abundance ratios between the protonated 
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ions and their neutral precursors. Then, a direct 
measurement of the protonated ion abundance leads 
to an indirect determination of the neutral abund- 
ance. A neutral abundance obtained by our analysis 
is ideally dependent on the following assumptions 
only: 

1) a cosmic ray induced ionization rate 

2) the occurrence of rapid exothermic ion- 
molecule reactions at interstellar cloud temperatures. 
However, the neutral densities we obtain are partially 
model dependent because certain results of the HK gas 
phase model (e.g., electron densities) must sometimes 
ba utilized. Note that the possiblity of reaction 

on grain surfaces has not been excluded in this paper. 

The motivation for this paper comes from the 
observation of HCO + by Buhl and Snyder (1970) , and 
Snyder et al. (1975) plus that of N 2 H + by Turner 
(1974), and Thaddeus and Turner (1975). These ob- 
servations demonstrate that polyatomic ions exist 
in appreciable abundance in the interstellar medium 
and that, consequently, ion-molecule reactions must 
occur. Until recently, identification of inter- 
stellar HCO~*' and N 2 H + , being based on theoretical 
arguments, (e.g.. Green et al. 1974) had been re- 
garded by some as equivocal. The recent success 
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of Woods et alo (1975) in observing the microwave 
spectra of these two ions in the laboratory has con- 
firmed their interstellar identification. In addi- 
tion,, the laboratory technique of Woods et al. (1975 
see also Saykally et al. 1976) promises to aid as- 
tronomers greatly in their search for a richer var- 
iety of interstellar molecular ions. 

JL 

Since N 2 H has now been unambiguously identi- 
fied in interstellar clouds, our method can be ap- 
plied quantitatively to the abundance of molecular 
nitrogen.. It is deduced that N 2 is probably the 
major nitrogen species in dense clouds. This point 
is noteworthy in view of the failure to observe N 2 
in diffuse clouds with the Copernicus satellite 
(Snow 1975) » For the other important neutral spe- 
cies, 0 2 ^ C 2 , C0 2 , H 2 0, C 2 ri 2 , CH^, observation of 
the protonated species has not yet been achieved in 
either the laboratory or interstellar clouds. In 
this paper estimates for the rotational transition 
frequencies of some of these ions are provided along 
with discussions of their probable column densities. 
It is to be hoped that highly accurate experimental 
frequencies of some of the ions discussed here will 
soon become available so that extensive searches for 
the protonated species of unobservable neutrals can 
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be undertaken.* 


II. DIATOMIC SPECIES 
a) N 2 


To determine the molecular nitrogen abundance 

+ . . 

from the observed N 2 H column density,, it is nec- 
essary to consider the ion-molecule chemistry re- 
lating to N 2 and N 2 H + „ Table 1 shows the relevant 
reactions together with their rate constants. It 
is assumed that CO is much more abundant than C0 2 
and fLjO. If these latter molecules are, on the 
contrary, as abundant or even more abundant than 
CO, more reactions must be considered in the analy- 
sis. However, the inferred N 2 abundance is sensi- 
tive only to the Sum of CO, C0 2 , and H 2 0 concen- 
trations, and is virtually independent of their 
relative abundances. Dissociative recombination 

•j- , t t t 

of N 2 H with electrons is ignored m our analysis. 
If the CO concentration reflects a strongly depeted 
carbon elemental abundance rather than the normal 
cosmic abundance ratio assumed in HK, and if the 
electron abundance is an order of magnitude higher 
than the HK values, then electron-ion recombina- 
tion cannot be. ignored. 
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4 - 

We cannot make the assumption that N 2 H comes 
from N 2 because of the likelihood thcit reactions 

4 * •}< 

such as + N = N 2 H + H 2 occur. However, we 

can obtain a lower limit to the molecular nitrogen 
density by assuming that the only mechanism for N 2 
formation involves N 2 H + and CO (both unambiguously 
observed species): N 2 II* i " + CO HCO**" + N 2 • Once 

formed, N 2 is almost' an inert species, reacting with 
no neutrals known to be present in abundance. Mo- 

♦ 4 4 - 

lecular nitrogen does react with the ions H 3 and 

He -1 *, formed via cosmic ray bombardment (Table X) 0 

1 

The steady-state assumption applied to N 2 leads 

i 

Assumption of the standard grain parameters leads 
to the conclusion that the ion-controlled depletion 
rate of N 2 is probably slower than the maximum con- 
densation rate on grains. If we assume that the ISL, 
absorbed on grains is unreactive, then the only ef- 
fect of such a depletion process is to transfer 
some of the gas phase N 2 onto the grains. 


to the relation (valid if [2ST 2 J <<: [CO]) 
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I5d. a h 

[CO] 2 k 4a k la /k 5a + 0,20 k 4b k lb /k 5b [H 2 i 

where the cosmic abundance of helium is set equal 
to 0.14 of the hydrogen abundance (Allen 1973). The symbol. 
"[ ]" stands for density and the k.j_ refer to the 
appropriate rate coefficients. Turner (1974) has 
estimated that i"- the clouds he observed [N 2 H ]/ 

[il 2 ] » 3 x 10 -11 ' with an order of magnitude uncer- 
tainty. For the CO density in relation (l)j it is 
reasonable to substitute the cosmic abundance of 
carbon multiplied by twice the H 9 density (HK) . 

Then., for a "typical” dense cloud in which [H 2 J ~ 

4 -3 

10 cm „ the inferred nitrogen density is [N 2 3 & 

-3 

0.1 cm . An upper limit to the N 2 density is given 
by cosmic; abundance considerations (HK) j viz.j [N 2 3 ^ 

A 

1 era similarly; for a cloud with [H 2 3 ~ 3 x 10 

_3 _3 ^3 

cm ^ 1 cm ^ [N 2 3 ^ 3 cm . It is doubtful that 
higher density regions were sampled by Turner (1974). 
Snyder (1975) has reported an observation of dim- 
xnishmg N 2 H density towards the denser regions of 
Orion Aj but detailed numerical figures are not yet 
available. The proximity of the lower and upper 
limits to the N 2 concentration is strong evidence 
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that, at the least, a significant fraction of the 

nitrogen cosmic abundance exists in the form of & 2 

4-3 

in clouds of density n £ 10 cm <> This conclusion 
also obtains in the full HK model. If a significant 
portion of the carbon were depleted from the gasi . 
phase,, the quantitative but not the qualitative con- 
clusions of this analysis would nave to be modified. 

It is interesting to speculate on the atomic 
vs. molecular fraction of the nitrogen abundance. 

The intensity of the observed N 2 H emission appar- 
ently suggests via the above analysis, that to 
within an order of magnitude all nitrogen is in the 
diatomic form. (The observed concentration of 
ammonia is too small to be of significance in this 
discussion. ) Yet, the uncertainties in both ob- 
servation and analysis certainly do not preclude 
the statement that atomic nitrogen may be of equal 
or even greater abundance than N 2 . Indeed., the 
existence of cosmic ray-produced helium ions will 
always serve to produce some N from N 2 (Farragher 
1970, Laudenslager et al. 1974) : n 2 + He + -* + n 

+ He. Thus, only in regions (if they exist) suffi- 
ciently dense to exclude 100 MeV cosmic rays is it 
feasible to neglect atomic nitrogen. 
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b) 0 2 


Several group* have studied the kinetics of 
the reaction H 3 + + l> 2 «- 0 2 H + + H 2 (Fehsenfeld, Lind- 
inger,. and Albritten 1975; Kim, Beard, and Huntress 1976). 
It appears that the reactants and products are iso- 
energetic, At 298°K, both the forward and reverse 
reactions are rapid with K, the equilibrium constant, 
near unity. ' . . in addition, the forward and 

reverse rate coefficients are virtually independent 
of energy. It therefore must be tentatively con- 
cluded that K si 1. at interstellar cloud tempera- 
tures, Because depletion by reaction with molecular 

hydrogen is so efficient, no other depletion mech- 

+ 

anisms need be considered for 0 2 H , There are no 
other important formation mechanisms than the one 
mentioned above. Therefore, at -steady state. 


[0 2 h + ] 


= [H 3 + ] 


[°2 ] 

Iv 


k l 

k 5 [CO] 


( 2 ) 


J r 

To obtain a severe upper limit to the 0 2 H density, 

assume that [0 2 ] [CO], Then [0 2 H + ] < 1.0 x 10 8 

-3 19 

cm . For a cloud of dimension 10 cm, we can 

+ . 11-2 
expect a maximum 0 2 H column density of 1.0 x 10 cm , 



v/aLl below the sensitivity limit of current obser 
vational technique. Thus, our method cannot yet 
bo utilised as an indirect means of 0 2 detection. 
In general, any ion that reacts rapidly with H 2 
v/ill have an abundance too low to be measured at 
this time. 


c) C 2 

The formation of C 2 in interstellar clouds 
of moderate density has been discussed by Solomon 
and Klemperer (1972) . The reaction of C 2 with 

X t *[* t , t 

H 3 to give C 2 H is exothermic and is assumed to 

-t- . 

be rapid. C 2 H can also be formed via the se- 
quence of reactions: 

C + + CH - C 2 + + H 

C 2 + + H 2 “* C 2 H+ + H * 

Unfortunately, analogous to the case of 0 2 dis- 
cussed above, cannot be used to monitor the 

C 2 abundance since it reacts rapidly with H 2 form 

t 4* * » « *f* 

mg C 2 H 2 j an< ^ e( 3 ull;L kr rum abundance of C 2 I-I 

is therefore expected to be too low to detect 
(Kim, Beard, and Huntress 1975) . In fact, this 

reaction sequence terminates with dissociative 

, . + . 
recombination of C 2 H 2 to give C 2 H, 
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a known interstellar molecule (Tucker, Kutner, 
and Thaddeus 1974) . 


III. TRIATOMIC SPECIES 
a) C0 2 


The molecular ion HOCO rs formed and depleted 
principally by those reactions sha^n in Table 2. 

4 - 

The steady-state abundance of HOCO is then 


r ( k 6 [H 3 J * k 7 [N 2 H 

J ~ l k 8 J [CO] * 


(3) 


Reaction (8) has not been investigated experimen- 
tally. If it were unexpectedly slow, or if CO 
were not the dominant non-hydrogenic neutral spe- 

4 - 

cies, other depletion mechanisms for HOCO might 
gain importance; viz., 

HOCO + + H 2 0 -* H 3 0 + + C0 2 

HOCO + +0 - HCO + + 0 2 

H0C0 + + N "* BIC0 + + NO 

HOCO' 7 * + e ^ C0 2 + Ho 

JL , 

If these HOCO reactions with neutrals occur 
as expected at the normal Langevin rate, inclusion 
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of species such as ^0, 0 and N in the analysis 
would only affect tha result to the extent that 
CO in Eq. (3) would have to be replaced by the 
sum of several densities. If [CO 2 ] an< ^ as 

expected, are « [GO], Eq. (3) simplifies to 

[CO,] k-k ft [HOCO’ 1 '] 7 [H0C0 + 3 

— « 7.7 x 1.0 ' — (4) 

[CO] la 6 [Hg] [H 2 ] 

Green et al. (1976) have recently calculated 
the electronic structure and equilibrium geometry 
of HOCO . According to their calculation, the 
structure is that of a slightly asymmetric prolate 
top, with the rotational constants listed in Table 
3; also listed in this table are the predicted 
frequencies of several a-type microwave transi 7 
tions betweem low-lying (k_ 1 = 0) rotational levels. 
These rotation constants and transition frequencies 
are reasonably insensitive to minor adjustments in 
the equilibrium geometry, and the uncertainties 
quoted in Table 3, which are slightly less than 1%, 
are hopefully conservative. 

The asymmetrical geometry of HOCO' has an 
interesting virtue from an astronomical standpoint: 
it endows the rotational spectrum with a specificity 
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not found in a linear molecule, and makes an astro- 
nomical search worthwhile even in the absence of 
precise line frequencies. Consider for example the 
J — 4-3 transition, which is of interest because 
it will probably be stronger in space than the 
other lines in Table 3, and because it happens to 
lie in a part of the spectrum where extremely sen- 
sitive receivers currently exist. In addition to 
the K_^ = 0 transition listed in Table 3, there 
are nearby in frequency three K-type doublets cor- 
responding to K_^ =1, 2, and 3. The upper levels 
of these transitions lie respectively 44, 144, and 
310 K above the rotational ground state, so these 
transitions will definitely be excited in compara- 
tively hot molecular clouds like Ori A. 

This rotational structure constitutes a highly 

• • , *{• 
distinctive fingerprint for HOCO . The = 1 

doublet is widely split, by 4(B Q - C q ) » 650 MHz j its 

center of gravity will fall near, but (because of 

centrifugal distortion) not precisely on the K_^ 

- 0 line. To judge from the centrifugal distortion 

parameters of the similar molecules HNCO and HN^, the 

centrifugal shift of the center of gravity should be 

about 6 MHz to lower frequency. The = 2 doublet, 

on the other hand, is expected to be barely resolved: 


1.6 


the calculated splitting is 0.28 MHz,, which is at 
the limit of resolution in a very narrow line source 
like OMC-2j the center of gravity of this doublet is 
expected to be centrifugally shifted about 24 MHz to 
lower frequency. Finally,, the k_^ — 3 doublet will 
be totally unresolved,, and is expected to be shifted 
about 54 MHz to lower frequency. It is clear that 
even if some of these lines could be detected enough 

information would be available to almost certainly 

. . + 
identify the HOCO molecule. 

With these considerations in mind, we have at- 
tempted to find HOCO* with the NRAO 36-foot telescope*. 


*The National Radio Astronomy Observatory is operated 
by Associated Universities, Inc., under contract to 
the National Science Foundation. 


Two of the best molecular sources, Ori A and DR21(OH), 
were systematically surveyed in frequency over a 
range which exceeded the uncertainty in the J = 4-3 
transition according to Table 3 (the equally good 
source Sgr B2 was unfortunately too near the sun at 
the time of this survey to be observed) . In particu- 
lar, the frequency region 83.25 to 87.05 GHz was sur- 
veyed in Ori A, and the region 80.40 to 86.875 was 
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surveyed in DR21(0H). No lines were detected in 
either source which could be attributed to HOCO . 

A uniform upper limit cannot be precisely set be- 
cause the sensitivity of the search was somewhat 
different at different frequencies^ but it is 
probably safe to conclude that any HOCO line 
stronger than 0.5K (radiation temperature) would 
have, been picked up in either source. 

This upper limit can in turn be converted into 
interesting upper limits on the C0 9 abundance and 
the CO 2 /CO ratio. It will be assumed that the 
observational limit refers to the I<_^ = 0 transi- 
tion i.e. j the transition ^o4’'^03- s an< ^ we con ~ 

fine our attention to Ori A. For an optically 
thin line of radiation intensity T the column 

JR 

density is 


Srrv* AvkT 


N HOCO + 


R 


hc^Af 


( 5 ) 


-5 -1 

where A = 1.31 x 10 sec is the Einstein co- 
efficient of the 4 q4 level calculated from the ab 
initio value of 2 D for the permanent electric di- 
pole moment along the 0C0 axis (Green et al. 1976) j, 
and f is the fractional population of the 4 q 4 level. 
A reasonable lower limit on f is 0.05j based on the 
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J, 

excitation of molecules similar to HOCO in Ori 
A. Then taking Av = 4 km sec"' 5 ' for the linewidth 

JL 

in Ori A (the observed linewidth of HCO ) , and T R 
< 0.5 K, we find from Eq. (4) that 

%oco + S- 4 x l ° 13 om " 2 - <s> 

+ 

The Ori A source in HCO and other ions is 
18 

probably ~ 10 cm in the line of sight, so for 
a mean HOCO density we find 

[HOCO + ] < 4 x 10~ 5 . (7) 

Finally, via Eq. (4) it is possible to set a 

crude upper limit on the C0 2 /C0 ratio in Ori A. 

Taking as reasonable values for this source [CO] ^ 
-3 23 -2 

10 cm , N h « 2 x 10 cm , and assuming that 
^ + 

the density ratio of HOCO to H 2 is equal to the 
ratio of corresponding column densities, we find 
that, 

[C0 2 ]/[C0] < 0.15 (8) 

Current notions of the formation of C0 2 in 
molecular clouds are too vague to warrant much 
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discussion of this limit. It is to be hoped that 
4. 

the HOCO frequencies will soon be refined by labora 
tory measurements so that more stringent astrophysi*- 
cal searches can be performed. 

b) H 2 0 

The recent observation of the isotopic spe- 
cies HDO (Turner et al. 1975) implies that H 2 0 is 
present in large abundance in dense clouds. Ob- 

JL 

servation of the H 3 0 (oxonium) ion would permit 

another * perhaps more precise, determination of 

the H 2 0 abundance. The oxonium ion is expected 

to be one of the most abundant molecular ions (HK) . 

It is ioselectronic with ammonia anv should possess 

a similar inversion-rotation spectrum. Novick, 

Stevens, and Klemperer (1975) Lave performed an 

* 4 * • » 

ab initio calculation of the H^O inversion po- 
tential. These authors estimate an inversion fre- 
quency of 84 cm“' L ; unlike the 24 GHz NH^ inversion 

+ 

frequency. If this calculation is correct, H 3 0 
posses no transitions favorable to microwave de- 
tection. 

The observation of HDO (Turner et al. 1975) 

7 —3 

in the Kleinmann-Low nebula ([H 2 ] ~ 10 cm ), 
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mentioned above, leads to speculation that gase- 
ous H 2 0 consumes much of the oxygen abundance 
in dense clouds. This possibility can probably 

J. 

be excluded for most HCO sources. Our reasoning 
is as follows: if H 2 0 were present in "appreciable" 

abundance (e.g., [H 2 0] [CO] ) , HCO would be de- 

pleted by the reaction, 

HC0 + + H 2 0 - H 3 0 + + CO 

more efficiently than by the sole HK mechanism, 

HC0 + + e - CO + H. 

The result would be that HCO + would lose its posi- 
tion as the dominant polyatomic ion and be reduced 
to a column density no larger than that of other 

ions such as N 2 H «, However , in almost all sources 
« 4* 

except Orion A, HCO is thought to' be .considerably 
more abundant than N. 2 H + (Turner 1975) t We there- 
fore conclude tentatively that in aid. livelihood, 

H 2 0 is not as abundant a species as CO in dense 
clouds. 


IV . POLYATOMIC SPECIES 
a) C 2 H 2 

Table IV contains the ion-molecule chemistry 

-f- 

of the vinyl cation C 2 H 3 * 


It is assumed that 
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protonation of acetylene (C 2 H 2 ) occurs predomi- 

•U 

nantly via reaction with HCO . The vinyl cation 
is not depleted by reaction with abundant CO, but 
is probably depleted by reaction with such species 
as 0, N, 0 2 , and H 2 0 as well as by electron recoin- 

•j- 

bination. The assumption that C 2 H 2 is at steady- 
state leads to the relation 


l C 2 H 2> = { 


k io teJ 


+ k 
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[h 2 o] 


+ k^ 2 [0] + ki^ [K’J + 



f C 2 H 3 + ] 

[HC0 + ] 


(9) 


The acetylene abundance can therefore be estimated 
from knowledge of the vinyl cation to HCO abun- 
dance ratio. To perform this estimation, we assume 
that, 

[H 2 0] + [0] + [Nj + [0 2 J « 2 x 1C)" 4 p 2 ]. 

This assumption is based on HK estimates but can 
be defended on simple cosmic abundance arguments. 
The dissociative electron recombination contribu- 
tion to CjHg ** depletion can probably be ignored 
in a rough analysis unless electron abundances 


v 


are somewhat higher than currently estimated. 

With these simplifications, relation (9) reduces 

> 

to 


[ c 2 h 2 ] 

[H 2 J 




2 x 10“ 4 


[HCO + ] 


( 10 ) 


4 * 

HCO 


column densities are not known because of 


.14-15 -2 


unknown opacities bu t N^ C Q + & 10 cm 

is a reasonable estimate (Turner 1975). A measure-* 


ment of N c H + » 10‘ L3 cm 2 would then imply that 

2 3 r— 

[C 2 H 2 I/[H 2 ] « 2 x 10" D - 2 x 10 . 


It should be noted that the vinyl cation cannot be 

4 . 

produced from C 2 H because C 2 H 2 > the protor\ated 

JL 

ion of C 2 H, does not react with H 2 to form CgHg 
(Kim, Beard, and Huntress 1975) . Therefore, 
is an indicator of acetylene density only. 

The structure of the vinyl cation has been 
studied’ theoretically by several groups (Hariharan, 
Lathan, and Pople 1972; Zurawski, Ahlrichs, and 
Kutzelnigg 1973; and Weber et al. 1976) 1 There 
are two possible structures of similar energy: 



and it is unclear which geometry possesses the 
greater stability, although the most reliable 


calculations seem to favor structure (b) . Struc- 
tures (a) and (b) are both slightly asymmetric 
prolate tops. Structure (b) permits only strong 
b-type transitions at frequencies in the far infra- 
red, whereas structure (a) allows AK = 0 and K- 
doubling a-type transitions. recent calculation 
{Weber et al. 1976) indicates that the barrier 
against tunneling between the two structures is only 
500 - 1500°K. If this range of values is correct, 
such a motion would have to be considered in an 

i 

analysis of the expected spectra, if interstellar 
detection of the vinyl cation is to be accomplished 

with current microwave receiver s_, structure {a) 
will have to be thermally populated at low tempera- 
ture. Table V lists the rotational constants for 
(a) as well as several rotational transition fre- 
quencies . Error estimates are uncertain but are 
at least +1 GHz for the Aj = 1 transitions. 

b) ch 4 

The abundance of methane in dense interstellar 
clouds is of great interest. Is gas phase carbon 
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mostly in oxidized species such as CO and C0 2 
or is a high percentage tied up in saturated hydro- 
carbons such as CH^? Observation of CH^* might 
help to elucidate this problem. There exist a 
number of theoretical studies of the structure of 
this ion (Lathan, Hehre, and Pople 1971' ; Hari- 
haran, Lathan, and Pople 1972; Dyczmons, Staemmler, 
and Kutzelnigg 1970) . Several structures — of C 4v , 
D 3h an< ^ C s s yn™e fcr y‘-have been investigated. Of 
these, the asymmetric C s structure appears to be 
the lowest in energy by less than 1 eV. Guest , 
Murrell, and Pedley (1971) found the potential 
surface of CH^ to be extremely flat at ;the 
equilibrium geometry. These authors also deter- 
mined that two of the hydrogens could undergo 
virtually free rotation. The expected microwave 
spectrum of such a poorly defined structuce with 
large amplitude motion is difficult to deduce. 

Much more theoretical work is needed. 


V. OTHER IONS OF IMPORTANCE 

In the preceeding section, the observation 
of protonated ionic species has been discussed 
as a means of obtaining information about unobservable 
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neutrals. There are other ions which are of 
interest to an understanding of interstellar cloud 
chemistry for a variety of reasons. Two of these 
species merit special concern. 

The ion H 2 CN is thought to be a precursor 
of IiCN (HK, Watson 1974, H&rbst 1976) . There is 
reason to believe that the HK model underesti- 
mated the abundance of this ion. Pearson and 

Schaefer (1974) have performed ab inition calcu- 

+ 

lations on HCNH and have deduced a linear struc- 
ture with B q = 37.4 GHz. The uncertainty in the 
prediction of the J = 1 -* 0 transition, lying at 
74.8 GHz, is probably less than 1 GHz, and may be as 
small as 300 MHz. 

Herbst (1976) has suggested that the ion 

4. 

H 2 CNH 2 can be formed by an especially favorable 
radiative association process, 

h 2 cn + + h 2 - h 2 cnh 2 + + hv . 

4. 

Methylenimmonium (H 2 CNH 2 ) can then recombine 
dissociatively to form H 2 CNH (methanimine) or 
add another H 2 via radiative association. Observa- 
tion of H 2 CNH 9 (especially coupled with observa- 

4. 

tion of HCMH ) could test this hypothesis. Kollman 
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4 * 

et al. (1973) have calculated a structure for HjCNH . 
(See Table VI.) Based on this structure, we have 
computed the rotational constants of this slightly 
asymmetric prolate top and tabluated several ro- 
tational transition frequencies of interest (Table 

4 * 

VI) . As is the case for C 2 H 3 > t ^ e uncertainty 
in the frequencies of the AJ = 1 lines is likely 
to be at least + 1 GHz . 

VI. DISCUSSION 

In this paper we have considered an indirect 
method for the determination of molecules that are 
normally unobservable under interstellar conditions. 
The ion-molecule kinetics and rotational spectra of 
the protonated ions of N 2 .> 0 2 , C0 2 , H 2°' t C 2 H 2 anc ^ 

CH^ have been discussed. Since laboratory microwave 
absorption spectra of ions has been demonstrated by 
Wood^ our estimates of transition frequencies hope- 
fully will be replaced by measurements. Of consid- 
erable importance in understanding the molecular 
composition of the interstellar medium is the active 
oxygen to carbon ratio. If indeed most molecular 
formation occurs in the gas phase, then the 0/C ratio 
in that phase should be reflected by the molecular 
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composition observed by radio astronomy. In this 

article we have suggested that present observations 

show both H 2 O and C0 2 are much less abundant than CO. 

We note that C0 2 is the expected form of excess oxygen., 

from equilibrium considerations. In particular the 

reaction C0 2 + H 2 ~ CO 4- H 2 0 is endothermic by 0.5 ev. 

The conclusion that C0 2 and H 2 0 are considerably 

less abundant than CO suggests that the gas component 

of the interstellar medium is more strongly depleted 

with respect to oxygen than carbon. This conclusion 

is based primarily on observations in dense clouds. 

The Copernicus observations, which are direct, have 

sufficient uncertainty that it is not possible to es~ 

tablish the 0/C gas phase ratio usefully in diffuse 

clouds. The high fraction of atomic oxygen (0 ) and 

X XX 

atomic carbon (C and C ) in diffuse clouds implies 
that the molecular abundances therein will not be 
seriously affected by the 0/C ratio.. 

The presence of molecular nitrogen is established 
by the abundance of N 2 H . The important question of the 
relative ratio of N 2 /N cannot be well answered by N 2 H 
alone. 

We have received a manuscript of a very rele- 
vant article by Snyder, Hollis and Watson, "An Inter- 

4* 

pretation of the Anomalous Variation of N 2 H /l-IC0/S0 2 
in Orion: Support for Ion-Molecule Reactions". 


This 
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article establishes limits on the abundance of H 2 0 
in agreement with the results discussed here. 
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TABLE I 

Major Ion-Molecule Reactions 

Reaction 

1 a) H 2 + C.R. Proton - H 2 + + e + C.R.P.' 
b) He + C.R. Proton -• He + + e + C.R.P.' 

2 H 2 + + H 2 - H 3 + + H 

3 N 2 H + + CO - HCO + + N, 

4 a) N 2 + H 3 + - N 2 H + + H 2 

b) II 2 + He + -• Products 

5 a) CO + H 3 + - HCO + + 

b) CO + He + - C + + 0 + He 


Involving N-,H + and N 2 


Rate Coefficient (k) 

Reference 

0.95 x 10~ 17 S _1 

Herb3t & 

1.0 x 10~ 17 S _1 

Klemperer 1973 

2.1 x 10 _9 cm 3 S" 1 

Bowers et al. 
1969 

8.8 x 10 -I0 cm 3 s“ 1 

Herbst et al. 
1975 

1.8 x 10" 9 cm 3 S _1 

Bohme et al. 
1973 

-9 3 -1 

1.7 x 10 crnS 

Kim et al. 
1976 

-9 3 —1 

1.2 x 10 cinS 

Farragher 

1970 

Laudens lager 
et al. 1974 

-9 3-1 

1.4 x 10 cm S 

, -9 3-1 

Burt et al. 
1970 

1.7 x 10 cm S 

Kim et al. 
1976 

1.7 x 10* 9 cm 3 S _1 

Fehsenfe Id 
et al. 1966 

-9 3 -1 

1.6 x 10 cm S 

Laudens lager 
et al. 1974 


keproduci v ;-uty of tub 
t>RKTXKAT : . -S 18 POOR 
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TABLE II 


Major Ion-Molecule Reactions Concerning HOCO 


+ 


Reaction Rate Coefficient 

6 C0 2 + H 3 + - HOCO + + H 2 1.9 X 10“ 9 cm 3 S" 

4* j. —10 3 

7 C0 2 + N 2 H - HOCO + N 2 9.2 x 10 cm S 

8 HOCO + + CO - HCO + + C0 2 1. X 10 -9 cm 3 S" 1 


Reference 
Burt et al. 

Burt et al. 


1970 
19 7 O' 


Estimate 
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TABLE III 
HOCO* 

A) Rotational Constants (MHz) [good to better than 1%] 

A = 706,000 
o 

B q = 10,820 

■C = 10,660 

b = - 1.15 x 10" 4 

JP 

B) Rotational Transition Frequencies 

J' K' _ 1 J" K"_ 1 K" x u (MHz) 

1 0 
2 0 
3 0 


1 000 21,480 + 200 

2 101 42,960 + 400 

3 202 64,440 + 600 


4 0 


4 


3 0 


3 


85,920 + 800 


A. ■ 
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table IV 


Major Ion-Molecule Reactions Concerning C 2 H 3 


+ 



Reaction 


Rate 

Coefficient 

9 

C 2 H 2 + HC 0 + A C 2 H 3 + + CO 

2 

X 

10 -9 cm 3 

s- 1 

10 

C 2 H 3 + + 

e -* Products 


2 

X 

10~ 6 cm 3 

s -1 

11 

c 2 h 3 + + 

H 2° ^ C 2 H 2 + 

4* 

h 3 o 

1. 

11 

-9 3 

X 10 . cm S 

12 

C 2 H 3 + + 

0 - H 2 C 2 0 + + 

H 

2 

X 

10 _9 cm 3 

s' 1 

13 

c 2 h 3 + + 

N - h 2 C 2 N+ + 

H 

2 

X 

10 ~ 9 cm 3 

s- 1 

14 

C 2 H 3 + + 

°2 - C 2 H 2°2 + 

+ H 

2 

X 

_q 3 
10 J cm J 

s" 1 


Reference 

Estimate 

Herbs t and 
Klemperer 1973 

Huntress ( 1976 ) 

Estimate 

Estimate 


Estimate 
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TABLE V 


H 2 CCH + (Vinyl Cation) 


A) Rotational Constants (MHz) 

A Q = 282,000 

B = 33,250 

. o 

C = 29,750 

o 

b = -7 X 10" 3 

.P 


B) Rotational Transition Frequencies 


J' K'^ 

10 1 

2 0 2 

110 
2 11 


J" K , '_ 1 K’ ^ 

0 0 0 
10 1 

2 12 


u (MHz) 

63,000 

126,000 

3,500 

10,500 
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TABLE VI 

H 2 GNH 2 (Methylenimmonium) 


A) Rotational Constants (MHz) 


A Q « 152 , 000 

B = 29,900 

o 

C « 25,000 

o 

b p - - .02 


B) Rotational Transition Frequencies 



u (MHz) 

54,900 

109,700 

4,900 

14,700 
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Pig. 1 


Pig. 2 


FIGURE CAPTIONS 


+ 

Calculated Structure of HOCO from Green 
et al. (1976) . 

Structure of H 2 C 2 H (vinyl cation) computed 
by Hariharan et al. (1972) . 

Structure of fLjCNEt, (Methylenimmonium) 
computed by Kollman et al. (1973). 
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